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ABSTRACT 

We discuss a method to determine orbital properties and masses of low-mass bodies orbiting eclipsing binaries. 
The analysis combines long-term eclipse timing modulations (light-travel time or LTT effect) with short-term, 
high-accuracy astrometry. As an illustration of the method, the results of a comprehensive study of Hipparcos 
astrometry and over a hundred years of eclipse timings of the Algol-type eclipsing binary R Canis Majoris are 
presented. A simultaneous solution of the astrometry and the LTTs yields an orbital period of P12 = 92.8 ± L3 yr, 
an LTT semiamplitude of 2574 ± 57 s, an angular semi-major axis of 012 = 117 ± 5 mas, and values of the orbital 
eccentricity and inclination of ei2 = 0.49 ±0.05, and /12 = 91.7 ±4.7 deg, respectively. Adopting the total mass of 
R CMa of M12 = 1 .24 ± 0.05 Mq, the mass of the third body is M3 = 0.34 ±0.02 Mq and the semi-major axis of its 
orbit is fl3 = 18.7 ± 1 .7 AU. From its mass, the third body is either a dM3-4 star or, more unlikely, a white dwarf. 
With the upcoming microarcsecond-level astrometric missions, the technique that we discuss can be successfully 
applied to detect and characterize long-period planetary-size objects and brown dwarfs around eclipsing binaries. 
Possibilities for extending the method to pulsating variables or stars with transiting planets are briefly addressed. 

Subject headings: stars: individual (R CMa) — astrometry — stars: fundamental parameters — binaries: 
eclipsing — stars: late-type 



1. INTRODUCTION 

Within the next decade several space astrometry missions, 
FAME and DIVA then SIM and GAIA, capable of sub-milli- 
arcsecond to micro-arcsecond accuracy are expected to be 
launched. One of the primary scientific goals of these mis- 
sions is the astrometric detection of low mass objects around 
nearby stars, including brown dwarfs and Jupiter-sized planets. 
The detection of these objects will be accomplished through 
the observation of the reflex motion of the host star caused by 
' the gravitational pull of the low-mass body. Although these 
, missions are capable of very high astrometric accuracies and 
' can observe up to millions of stars, their lifetimes are relatively 
short (2.5 to 5 yr). Thus, these space missions are optimized to 
, detect planets within the habitable zones of late-type stars, but 
' they could fail to detect (additional) planets with longer peri- 
ods. It is important, however, to secure a complete picture of 
the bodies orbiting a star both from a pure census point of view 
and also to understand the genesis and evolution of planetary 
systems. In addition, planets do not necessarily remain within 
the habitable zone because of long-term chaotic perturbations. 
As we know from our Solar System, the presence of massive 
planets, such as Jupiter and Saturn, in distant orbits play a cru- 
cial role in stabilizing the orbits of the inner planets. 

One effective way of extending the time baseline that per- 
mits the discovery of long period exosolar planets or brown 
dwarfs is to use the light-travel time (hereafter LTT) effect in 
eclipsing binaries. From this technique, the eclipses act as an 
accurate clock for detecting subtle variations in the distance to 
the object (this is analogous to the method used for discover- 
ing earth-sized objects around pulsars, see Wolszczan & Frail 
1992). The periodic quasi-sinusoidal variations of the eclipse 
arrival times have a very simple and direct physical meaning: 



the total path that the light has to travel varies periodically as 
the eclipsing pair moves around the barycenter of the triple sys- 
tem. The amplitude of the variation is proportional both to the 
mass and to the period of the third body, as well as to the sine 
of the orbital inclination. As discussed by Demircan (2000), 
nearly 60 eclipsing binaries show evidence for nearby, unseen 
tertiary components using LTT effects. A recent example of a 
brown dwarf detected around the eclipsing binary V47 1 Tau us- 
ing this method was presented by Guinan & Ribas (2001). Also 
this method is being employed in selected low mass eclipsing 
binaries to search for extrasolar planets (Deeg et al. 2000). 

The primary advantages of using the LTT effect to detect 
third bodies in eclipsing binaries are: /) The necessary pho- 
tometry can be secured with small telescopes using photoelec- 
tric or CCD detectors; //) The number of eclipsing binaries is 
large - 4000 currently known in the Galaxy - and this number 
could increase very significantly when results from upcoming 
astrometry and photometry (e.g., MONS, COROT, Eddington, 
Kepler) missions are available; For select eclipsing bina- 
ries (with sharp and deep eclipses) the timings can be deter- 
mined with accuracies as good as several seconds; iv) The mass 
of the eclipsing pair can be known from conventional spectro- 
scopic and light curve analyses. A shortcoming of the LTT 
method is that only upper limits to the mass and size of the 
orbit of the tertiary component can be determined (the analysis 
yields the mass function"*, /(M3), and as sin /3). However, as it 
was demonstrated in the case of Algol (Bachmann & Hershey 
1975), the LTT analysis can be complemented with astrome- 
try to yield the orbital inclination and thus the actual mass and 
semi-major axis of the third body. Furthermore, with the or- 
bital properties (P, e and uj) known from the LTT analysis, only 
a small fraction of the astrometric orbit needs to be covered 
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when using high-accuracy astrometry. 

In this paper we present the results of the combined LTT 
analysis and Hipparcos astrometry of the Algol-type eclipsing 
binary R CMa. The residuals of over 150 eclipse timings ex- 
tending from 1887 to 2001 show a periodic (^93 yr) quasi- 
sinusoidal modulation. As previously shown by Radhakrish- 
nan, Sarma, & Abhyankar (1984) and Demircan (2000), these 
variations are best explained by the LTT effect arising from the 
gravitational influence of a third body. The Hipparcos astrome- 
try also shows the presence of small but significant acceleration 
terms in the proper motion components explicable by the reflex 
motion from a third body. Our study illustrates that with a well- 
defined LTT effect, only a few years of accurate astrometry are 
needed to constrain the orbital solution and determine the mass 
of the third body. 

2. OVERVIEW OF R CMA 

R Canis Majoris (HD 57167, HR 2788, HIP 35487) is a 
bright (Vinax = 5.67 mag), semi-detached eclipsing binary hav- 
ing an orbital period of 1.1359 days. As pointed out by Var- 
ricatt & Ashok (1999), R CMa holds special status among Al- 
gol systems in that it is the system with the lowest known total 
mass and hosting the least massive secondary star. Since the 
discovery of its variability in 1887 by Sawyer (1887), R CMa 
has been frequently observed and has well determined orbital 
and physical properties. The major breakthrough in under- 
standing the system came when Tomkin (1985) was able to 
measure the very weak absorption lines of the faint secondary 
star and determine the masses of the two stars from a double 
line radial velocity study. The analyses of its light and ra- 
dial velocity curves (see Varricatt & Ashok 1999) show that 
this system has a circular orbit and consists of a nearly spher- 
ical FO-1 V star (Mj = 1.07±0.2 M©; Ri = 1.48±0.10 Rq; 
Li/Lq = 5.78 ±0.38) and a low mass, tidally-distorted K2-3 
IV star (M2 = 0.17 ±0.02 M©; R2 = 1.06 ±0.07 R©; L2/L0 = 
0.43 ±0.10). Moreover, nearly every photometric study in- 
dicates that the cooler star fills its inner Lagrangian surface. 
The relatively high space motions (S = 67 km s~') suggest that 
R CMa is a member of the old disk population and thus a fairly 
old (5-7 Gyr) star (Guinan & lanna 1983). 

The present state of the system is best explained as a low 
mass Algol system that has undergone mass exchange and ex- 
tensive mass loss. Asymmetries in its light curves and subtle 
spectroscopic anomalies indicate that mass exchange and loss 
are still continuing but at a much diminished rate compared to 
most Algol systems. The very low mass of the secondary star 
and old disk age indicate that R CMa is near the end of its life 
as an Algol system. As in the case of all Algol systems, the sec- 
ondary star lies well above the main-sequence. However, unlike 
most Algol systems, the primary star is too hot and overlumi- 
nous for observed mass. Moreover, a recent analysis of older 
photometry of R CMa by Mkrtichian & Gamarova (2000) in- 
dicates that the F star is a low-amplitude S Scuti variable with 
a B light amplitude of 9 millimagnitudes and a period of 68 
minutes. 

3. OBSERVATIONS 

3.1. Astrometry 

Hipparcos observed R CMa between March 9, 1990 and 
March 5, 1993. There are 68 one-dimensional astrometric mea- 
surements corresponding to 35 different epochs in the Hippar- 
cos Intermediate Astrometric Data, which were obtained by the 



two Hipparcos Data Reduction Consortia (33 measurements 
from FAST and 35 from NDAC). The astrometric data can 
be obtained from CD-ROM 5 of the Hipparcos catalog (ESA 
1997). Unfortunately, the timespan of the Hipparcos obser- 
vations is much smaller than the orbital period of the tertiary 
component, and this might eventually result in possible sys- 
tematic errors in the orbital elements. To further constrain the 
solution additional older ground-based positions must be used. 
Indeed, Tycho-2 proper motions were computed by combin- 
ing Tycho-2 positions and ground-based astrometric catalogs. 
For R CMa, 17 epoch positions of ground-based catalogs used 
for the Tycho-2 proper motion computation were kindly made 
available to us by Dr. Urban and are listed in Table |l[ These 
measurements span over one century and so the Tycho-2 proper 
motions can be understood as the combination of the true proper 
motion and of a large fraction of the orbital motion. Conse- 
quently, the Tycho-2 proper motion of R CMa cannot itself be 
used in our analysis and only the individual positions contain 
valuable orbital information. In contrast, a short-term proper 
motion determination, such as the one computed around 1980 
by Guinan & lanna (1983), reveals itself to be very useful in 
constraining the astrometric solution. 

In the course of the astrometric data reduction of the Hip- 
parcos data, a test was applied to all the (apparently) single 
stars in order to check whether their motion was significantly 
non-linear. Most likely, a significant curvature of the photocen- 
ter motion is an indication of a possible duplicity. As it turns 
out, R CMa is one of the 2622 "acceleration" solutions of the 
Double and Multiple Stars Annex of the Hipparcos Catalogue, 
which provides a hint for the presence of a third body, indepen- 
dently from the LTT effect. 

3.2. Photometry 

R CMa has a long baseline of eclipse timings that extend 
from the present back to 1887. Most of the early eclipse times 
were determined from visual estimates. Several period stud- 
ies have been carried out. Early studies of times of minimum 
light indicated a possible abrupt decrease in the period during 
1914-15 (see Dugan & Wright 1939; Wood 1946; Koch 1960; 
Guinan 1977). However, as more timings accumulated it be- 
came apparent that the long-term variations in the (O-C)'s of 
the system are periodic and thus best explained by the LTT ef- 
fect produced by the presence of a third body. The analysis of 
Radhakrishnan et al. (1984), with eclipse timings from 1887 to 
1982, and Demircan (2000), who includes timings up to 1998, 
make a strong case in support of the LTT scenario. 

Our photoelectric eclipse timing observations extend the time 
baseline up to early 2001. The observations were obtained with 
the Four College 0.8-m Automatic Photoelectric Telescope lo- 
cated in Southern Arizona during 1995/96 and 2000/01. Differ- 
ential photometry was carried out using uvby Stromgren filter 
sets. The mid-times of primary minimum and the (O-C)'s for 
these are given in Table |, along with the corresponding uncer- 
tainties. The (O-C)'s were computed using a refined ephemeris 
determined from the analysis in §^ (Eq. Our observations 
were combined with those compiled from the literature to yield 
a total of 158 eclipse timings from 1887 through 2001. The pri- 
mary eclipse observations obtained from the literature are also 
provided in Table ^ Even though it is not explicitly mentioned 
in any of the publications, the times listed are commonly as- 
sumed to be in the UTC (coordinated universal time) scale. The 
timings were transformed from UTC to TT (terrestrial time) 
following the procedure described in Guinan & Ribas (2001), 
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which is based on the recommendations of Bastian (2000). The 
timings listed in Table ||are therefore HJD but in the TT scale. 

The uncertainties of the individual timings are difficult to es- 
timate. The compilation of Radhakrishnan et al. (1984) (from 
which most of the timings in Table ^ come) does not provide 
timing errors but only a relative weighting factor related to the 
quality of the data and the observation technique. We therefore 
adopted an iterative scheme to determine the actual uncertain- 
ties by forcing the ^ of the (O-C) curve fit (in §^ to be equal 
to unity. This rather arbitrary scale factor determination is in- 
deed justified because it ensures that the fitting algorithm will 
yield realistic uncertainties for the orbital parameters of the sys- 
tem. The individual timing errors are included in Table ^ The 
uncertainties yielded by the iterative scheme are about 10-13 
minutes for photographic timings and some 2-4 minutes for 
photoelectric timings, in both cases reasonable figures given the 
characteristics of the two methods and the shape of the eclipse. 

The Hipparcos mission, in addition to observing accurate po- 
sitions of R CMa, obtained a total of 123 photometric measure- 
ments, which are present in the Hipparcos Epoch Photometric 
Data (ESA 1997). The phase coverage of the observations is not 
sufficient to determine an accurate primary eclipse timing using 
conventional methods. As an alternative, we adopted the physi- 
cal information available to fit the entire light curve and derive a 
phase offset. Also, since the observations span 3 yr and the sys- 
tem exhibits (O-C) variations, the photometric dataset was split 
into two subgroups about January 1 , 1 99 1 . Then, using the pho- 
tometric elements of Varricatt & Ashok (1999), we employed 
the Wilson-Devinney program (Wilson & Devinney 1971) to 
run fits to both light curves by leaving only a phase shift and 
a magnitude zero point as free parameters. The fits were very 
satisfactory and we derived eclipse timings for the mean epochs 
from the best-fitting phase offsets. The resulting two timings, 
with uncertainties of around 100 s, are included in Table ^. 

These two timings are significantly different as one would 
expect from the LTT secular change of period during the 3- 
year duration of the Hipparcos observations. More precisely, 
each date of observation could be approximated by T « To + 
PoE + gE^, with E being the variability cycle number So, the 
photometric "acceleration" term g is a measure of the depar- 
ture from a linear ephemeris during the observation window 
(the Hipparcos mission lifetime in this case). Since Hippar- 
cos did not provide minimum timings, we used an alterna- 
tive method to fit the equation above. The epoch measure- 
ments were folded using the reference epoch and period from 
ground-based studies. Then, the quadratic term was estimated 
by minimizing the distance between successive points of the 
folded light curve (string-length method). The fit yielded a 
value of g = 2.09 • 10"*^ d, which leads to a cumulative effect of 
gE^.^^ = 0.0193 ± 0.0006 d during the course of the Hipparcos 
observations. This rough estimation gives a significant accel- 
eration term which is of the same order as that resulting from 
the long-term LTT analysis (§H) and indicates that the (O-C) 
variations attributed here to the presence of a third body could, 
in principle, have been detected through the Hipparcos photo- 
metric analysis alone. 



4. ANALYSIS 

The expressions that describe the LTT effect as a function of 
the orbital properties were first provided by Irwin (1952). In 
short, the time delay or advance caused by the influence of a 



tertiary component can be expressed as: 
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where c is the speed of light, and ai2, in, en, ton, and vn 
are, respectively, the semi-major axis, the inclination, the ec- 
centricity, the argument of the periastron, and the true anomaly 
(function of time) of the orbit of the eclipsing pair around the 
barycenter. As is customary, the orbital inclination in is mea- 
sured relative to the plane of the sky. The naming convention 
adopted throughout this paper uses subscripts 12 and 3 for the 
orbital parameters of the eclipsing pair and the tertiary com- 
ponent, respectively, around the common barycenter Obvi- 
ously, most of the parameters for the eclipsing pair's and the 
third body's orbits are identical (such as period, eccentricity, 
and inclination), but we use the subscript 12 because the actual 
measurements are made strictly for the brighter component of 
the long period system. Finally, subscript EB refers to the close 
orbit of the eclipsing binary. 

The fit of Eq. (jl]) to the timing data would provide a good es- 
timation of the a number of orbital and physical parameters of 
the system (see e.g., Guinan & Ribas 2001). However, both the 
orbital semi-major axis and the mass of the third body would 
be affected by a factor sin / that cannot be determined from the 
LTT analysis alone. When the LTT analysis is combined with 
astrometric data, all orbital parameters (including / and even 
Vi) can be determined yielding a full description of the system. 
The availability of Hipparcos intermediate astrometry permits 
the fitting of the observations using an astrometric model not 
accounted for in the standard Hipparcos astrometric solution. 
In the particular case of R CMa we have considered an or- 
bital model that has been convolved with the astrometric mo- 
tion (parallax and proper motion). The orbital motion produces 
the following effect on the coordinates: 
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In addition, since the Hipparcos measurements are unidimen- 
sional, the variation of the measured abscissa v on a great circle 
is: 

dv dv dv 

Av = -(Aacos(5 + Ajc) + — (A5+ Ax) + — AC3 + 

aa cos 00 aW 



9v . c dv ^ 
-A^Q cos(5+- — A^J.s 

a/ia cos Oflg 



(4) 



where the astrometric components are a, S for the coordinates, 
Ha, PI'S for the proper motion and tO for the parallax. 

In addition to fitting the orbital and astrometric properties of 
the system, a timing zero point and a correction to the orbital 
period of the eclipsing pair (that could lead to a linear secu- 
lar increase or decrease of the (O-C)'s) were also considered. 
The initial values of the period and the reference epoch were 
adopted from Varricatt & Ashok (1999). 

The full set of observational equations includes those related 
to the timing residuals (Eq. (|l])) and those coming from the as- 
trometric measurements (Eqs. (^, (||), & (^)). All these equa- 
tions were combined together and the 14 unknown parameters 
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(5 for the astrometric components - a, 6, ^a, f^s, Ci -, 7 for the 
orbital elements - 012, ei2, Lon, in, P12, 7^™' ^12 one for the 
reference epoch - Toeb ^nd one for the period of the eclipsing 
system - Peb -) were recovered via a weighted least-squares fit 
as described in, e.g., Arenou (2001) or Halbwachs et al. (2000). 
Note that the weights of the individual observations were com- 
puted as the inverse of the observational uncertainties squared 
and multiplied by the corresponding correlation factors. The 
uncertainties adopted are given in Tables |l] and || and in the 
Hipparcos Catalogue CD-ROM 5. 

Due to the short timespan of the Hipparcos observations, a 
large uncertainty on the reflex semi major-axis would exist if 
the Hipparcos astrometric data were used alone. As a first step 
towards better constraining the solution we added the epoch 
proper motion from Guinan & lanna (1983) as an external ob- 
servation (with the two subsequent equations for both compo- 
nents of the proper motion). To do so, the appropriate equa- 
tion for the first derivatives of the orbital motion was used. As 
expected, the quality of the solution improved and yielded a 
semi major axis of an = 140 ± 16 mas and a tertiary mass of 
M3 = 0.42 ±0.05 M0. Yet, a closer inspection of the resid- 
uals revealed that this solution was not fully compatible with 
the ground-based epoch positions mentioned in § [3.1| , as a clear 
trend appeared in declination. For this reason we decided to 
include these positions also in the fit, together with the Guinan 
& lanna proper motion and the photometric (O-C) minimum 
times. In total, the least-squares fit had 262 equations for 14 
parameters to determine. A robust fit approach (McArthur, Jef- 
ferys, & McCartney 1994) was used due to the large dispersion 
of the ground-based astrometric measurements. The resulting 
goodness of the fit was 0.63 and graphical representations of the 
fits to the eclipse timing residuals and Hipparcos intermediate 
astrometry are shown in Figure [l|. The less-accurate ground- 
based astrometric positions (with standard errors of about 200 
mas on average) are not represented in the figure for the sake of 
clarity. 

The resulting best-fitting parameters together with their stan- 
dard deviations are listed in Table || The astrometric solution 
presented supersedes that of the Hipparcos Catalogue because 
it is based upon a sophisticated model that accounts for the or- 
bital motion and considers ground-based astrometry. Also, Ta- 
ble H includes the mass and semi-major axis of R CMa C that 
follow from the adoption of a total mass for the eclipsing sys- 
tem. Finally, our fit also yields new accurate ephemeris for the 
eclipsing pair: 

r(Min I) = HJED2430436.5807 + 1 .13594197 E (5) 

where all times are in the TT scale and the zero epoch refers 
to the geometric center of the R CMa orbit. Note that the ac- 
curacy of the new period we determine (see Table |^) is better 
than 9 milliseconds. We have considered in our analysis a linear 
ephemeris as that in Eq. (^). However, Algol systems have been 
observed to experience secular decreases of the orbital periods 
possibly due to non-conservative mass transfer and angular mo- 
mentum loss (see, e.g., Qian 2000a). To assess the significance 
of this effect on R CMa, we modified our fitting program by 
considering a quadratic term. The coefficient of this quadratic 
term was found to be (-2 . 1 ± 1 . 1 ) ■ 1 0~' which translates into 
a period decrease rate of dP/dt = (-6.9 ± 3.6) • lO"*^ d yr"' . This 
is a very slow rate compared to other Algol systems (see, e.g., 
Qian 2001) yet commensurate with the low activity level of R 
CMa, which is near the end of its mass-transfer stage. Because 



of the poor significance (below 2cr) of the period decrease rate 
derived from the analysis, we decided to neglect the quadratic 
term and adopt a linear ephemeris. It should be pointed out, 
however, that the astrometric and orbital parameters resulting 
from the fit with quadratic ephemeris are well within one sigma 
of those listed in Table |3| 

Interestingly, a closer inspection of Figure [l| reveals small ex- 
cursions of the data from the LTT fit. To investigate these, we 
computed the fit residuals which are shown in Figure ^ The 
presence of low-amplitude cyclic deviations seems quite obvi- 
ous in this plot. If these (O-C) timing oscillations were caused 
by the perturbation of a fourth body in a circular orbit (R CMa 
D), its orbital period would be about 45 yr, with a LTT semi- 
amplitude of 275 s, a minimum mass of 0.06 M0, and an orbital 
semi-major axis of about 14 AU. The orbit of the third body is 
highly eccentric so that it would be interior to that of R CMa 
C near its periastron (r3p = 9.5 AU). The intersections of the 
two orbits would result in an apparently unstable configuration. 
Other possible explanations for the low-amplitude oscillations 
include abrupt period changes of the binary itself caused by 
variable angular momentum loss and magnetic coupling (see, 
e.g., Qian 2000b), a magnetic activity cycle of the secondary 
star (see, e.g., Applegate 1992), or simply a spurious effect 
caused by the inhomogeneity of the data set. Unfortunately, 
the available astrometric data are not sensitive enough to prove 
or refute the existence of a fourth body and only new accurate 
photometric eclipse timing determinations or high-accuracy as- 
trometry will provide the necessary evidence. 

5. DISCUSSION 

The orbital, astrometric, and physical properties presented in 
Table || are within 1 cr of the (less accurate) earlier estimates 
of Radhakrishnan et al. (1984), who based their analysis on 
eclipse timings up to 1982. However, our study, in addition to 
extending the time baseline, has been able to determine the in- 
clination of the third body's orbit by making use of the available 
high-precision astrometry (Hipparcos). Thus, R CMa joins Al- 
gol, the prototype of its class, in having the orbital properties 
of the third body determined from a combined analysis of the 
astrometry and LTT. The long period (^^93 yr) of R CMa C is 
the longest period detected and confirmed so far for an eclips- 
ing binary. This is chiefly because of the large LTT present in 
R CMa (total amplitude of 86 min) and the existence of eclipse 
timings available for this star back to 1887. 

It is interesting to note that the inclination of R CMa C is 
found to be of «92 ±5° and thus compatible with an edge-on 
value of 90°. In this situation, mutual eclipses of the tertiary 
component and the close binary pair might occur. This tantaliz- 
ing possibility is, however, very unlikely since eclipses are only 
possible within a very narrow window (^1.2 arcmin) about an 
inclination of 90°. If this were indeed the case, the transit of 
the tertiary component in front of the eclipsing pair should have 
occurred during mid 2001. Also interesting to note is the near 
coplanarity of the eclipsing system and its companion. Varri- 
catt & Ashok (1999) found an inclination for the eclipsing pair 
of /eb = 79.°5, which is equivalent to /eb = 100.°5 because of the 
degeneracy. Thus, the third body's orbit appears to be within 
only 8-13° of the orbit of the eclipsing pair. 

One question remains yet unaddressed, and this is the nature 
of the tertiary companion of R CMa. With a measured mass 
of 0.34 M0, one is tempted to classify R CMa C tentatively as 
a main sequence M3-4 star (Delfosse et al. 2000). However, 
another attractive possible scenario is a white dwarf (WD) as 
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tertiary component. There is no direct evidence for a WD com- 
panion to R CMa, but the mass of the third body is compatible 
with the low-end of the WD mass distribution found by Sil- 
vestri et al. (2001). The presence of a hot WD (Feft > 10000 K) 
is unlikely from lUE observations of R CMa in the UV region, 
where no hot source has been detected. Nonetheless, R CMa 
is an old disk population star so that a young WD is, in princi- 
ple, not expected. If a WD is present, its original stellar mass 
would have to be greater than the mass of the initial primary 
(now secondary) of R CMa. From binary evolution theory, the 
best estimate of the initial mass of the original primary is about 
1.4 Mq (Sarma, Vivekananda Rao, & Abhyankar 1996). This 
indicates a pre-WD evolution time for the companion of around 
2-3 Gyr Cooling sequences for WDs (Serenelli et al. 2001) 
yield an effective temperature of ^--^5400 K at an age of ^3 Gyr, 
which is a reasonable estimate given the kinematic characteris- 
tics of R CMa. Should the tertiary component turn out to be a 
WD, such an old and low-mass object might be exceedingly in- 
teresting since it could belong to the controversial class of blue 
WDs that have been claimed to play an important role in ex- 
plaining the dark matter content of the galactic halo (Hodgkin 
et al. 2000). 

As the predicted temperatures both in the WD and M-star 
scenarios are fairly similar, only the very different expected 
luminosities can help identify the nature of the companion to 
R CMa. Thus, the measure of the magnitudes of the tertiary 
component through direct imaging would be a definitive proof. 
If we consider the M-star scenario, the absolute magnitude of 
the tertiary component would be My ~ 11 mag, which trans- 
lates to my ~ 14 mag when using the parallax obtained in 
§Q This is about 8 mag fainter than R CMa itself. To give 
an example in the IR, the situation is significantly improved 
in the K-band, where the magnitude difference is reduced to 
AK w 4 mag. The tertiary component would be even fainter 
in the WD scenario. Indeed, the absolute magnitude can be 
estimated as My « 14 mag, which implies an apparent magni- 
tude of OTy « 17 mag. The difference with R CMa is therefore 
Ay « 1 1 mag. In the IR the situation does not improve signifi- 
cantly, with a large magnitude difference of AK « 10 mag. 

These magnitudes and dynamic brightness ranges are chal- 
lenging but yet attainable with state-of-the-art coronographs or 
Speckle spectrographs. Further complications arise from the 
current spatial location of the tertiary component near the con- 
junction of its orbit with the eclipsing pair Figure |^ depicts the 
predicted orbits of both the eclipsing pair and the tertiary com- 
ponent on the plane of the sky. As can be seen, the separation 
between the eclipsing system and R CMa C is only 27 mas as 
of 2002, which makes direct imaging very difficult. On an opti- 
mistic note, the situation will slowly improve in the future until 
a maximum separation of wO."8 is reached around year 2037. 

Claims of third body detections through the analysis of (O- 
C) residuals have sometimes been challenged. Spurious period 
changes caused by magnetic activity cycles, variable angular 
momentum loss, magnetic coupling, or other effects have been 
argued to explain modulations in the (O-C) residuals found in 
a number of eclipsing binary stars. Interestingly, R CMa would 
be a prime candidate for such spurious period changes because 
of its interactive nature. However, with over one period cy- 
cle in the LTT curve currently covered and, more importantly, 
with direct evidence from Hipparcos astrometry, the case for a 
tertiary companion to R CMa is now iron-clad. What only re- 
mains to be clarified at this point is whether this third star is 
an M dwarf or a WD. Also, the nature of the lower amplitude 



~45 yr variation needs to be further explored with continued 
observations. 

6. CONCLUSIONS 

This paper presents a combined analysis of short-term accu- 
rate astrometry and long-term timing residuals applied to the 
eclipsing binary R CMa. The study yields the complete orbital 
and physical properties of the tertiary component. A determina- 
tion of the mass of the third body is possible because the masses 
of the eclipsing binary components themselves are well-known 
from light and radial velocity curve analyses. 

The example discussed here illustrated the capabilities of a 
method that will reach its full potential with the upcoming high- 
accuracy astrometric missions. The improvements in precision 
of the future astrometric measurements are due to an increase 
up to a thousand-fold relative to Hipparcos and the quality of 
the photometry (and thus the eclipse timings) will also improve. 
More quantitatively, timings with accuracies of ~ 10 s are now 
possible for select eclipsing binaries with sharp eclipses. The 
detection of large planets (^10 Mj) in long-period orbits ('^10- 
20 yr) around eclipsing binaries will be therefore a relatively 
easy task. The short-term astrometry will confirm the detec- 
tions and yield the complete orbital solution (most significantly 
the inclination) and thus the actual mass of the orbiting body. 

One of the unexpected outcomes of the Hipparcos mission 
has been that a primarily astrometric satellite can also provide 
valuable new results from its photometric measurements alone 
(numerous new variables, HD 209458 planetary transits, etc). 
The data analysis of the next generations of astrometric satel- 
lites will surely benefit from a simultaneous analysis of the as- 
trometric and photometric data. Astrometric missions such as 
GAIA will likely detect one million new eclipsing binaries (a 
smaller number is expected for FAME). About one per cent of 
the eclipsing binaries observed by Hipparcos has a 0.0001 day 
precision in the reference epoch, which is enough to detect the 
LTT effect that would arise from a 10 Jupiter mass third body 
with all year period. If we assume the same ratio for GAIA, 
hundreds to thousands of third bodies would be detected. Al- 
though GAIA astrometry alone will be able to give the orbit for 
the closest stars, the orbit for more distant stars will depend on 
the availability of ground-based light curves to define the refer- 
ence epoch. 

This method of combining LTT analysis and astrometry com- 
plements very well with the ongoing spectroscopic searches. 
The LTT analysis favors the detection of long-period third bod- 
ies around eclipsing binaries because the amplitude of the time 
delay due to the LTT effect is proportional to P^^ while the 

-1/3 

spectroscopic semi-amplitude is proportional to Pjj . When 
the samples of spectroscopic and LTT systems are sufficiently 
large, we will have a complete picture of the distribution of 
bodies in a stellar system and a realistic test of planet formation 
theories will be possible. 

Finally, the LTT analysis method does not have to be nec- 
essarily applied to eclipsing binaries. In essence, the method 
is based upon having a "beacon in orbit", which, in the case 
of eclipsing binaries, are the mid-eclipse times. However, any 
strictly periodic event that can be predicted with good accuracy 
could be potentially useful to detect stellar or sub-stellar com- 
panions. This includes, for example, pulsating stars. More in- 
terestingly, transiting planets are also prime candidates for LTT 
studies. In this case, not only further orbiting planets could be 
discovered, but also good chances for detecting moons around 
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the transiting planet exist. 
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TABLE 1 

Ground-based astrometric data for R CMa. 

References for Table |l]: WASH - Urban (2001; priv. comm.); A.C. - Urban et al. (2000); CAPE - Urban (2001; priv. comm.); 
ALBANY - Urban (2001; priv. comm.); YALE - Urban (2001; priv. comm.); CPC2 - Zacharias et al. (1992); PERTH - Urban 
(2001; priv. comm.); CAMC - Fabricius (1993); FOKAT - Bystrov et al. (1991). 



Julian year a (deg) 5 (deg) ctq cos 5 (mas) ag (mas) Source 



1894.7 


109.86244875 


-16.391100556 


861 


690 


WASH AG 1900 


1905.15 


109.8630725 


-16.391588056 


253 


269 


A.C. 


1914.11 


109.863289583 


-16.391985556 


132 


116 


CAPE 2ND FUND 1900 


1914.97 


109.863497083 


-16.391931111 


304 


332 


A.C. 


1916.2 


109.863377917 


-16.392131667 


294 


282 


ALBANY 10 


1918.5 


109.863659583 


-16.392145833 


321 


270 


WASH 20 


1923.05 


109.863762083 


-16.392368333 


156 


153 


CAPE 1-25 


1933.52 


109.864227083 


-16.392731944 


316 


264 


YALE 12/1 -14/-18 


1934.44 


109.86428875 


-16.392647222 


237 


210 


CAPE 3-25 


1939.58 


109.86449 


-16.392880556 


130 


132 


CAPE 1-50 


1942.02 


109.864635833 


-16.393026944 


151 


200 


WASH 40 9-IN 


1969.33 


109.865958333 


-16.394058889 


45 


45 


CPC2 


1983.05 


109.866597083 


-16.394559722 


168 


164 


WASH TAC 


1984.37 


109.866713333 


-16.394581389 


86 


181 


PERTH 83 


1986 


109.8667525 


-16.394728333 


88 


119 


CAMC Series 


1986.07 


109.866875417 


-16.394721944 


110 


119 


FOKAT 


1992.49 


109.8670625 


-16.394969444 


39 


46 


WASH 2-JOO 
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Table 2 

Primary eclipse timings for R CMa. 

References for Table|: (1) Radhakrishnan et al. (1984); (2) Wood (1946); (3) Knipe (1963); (4) Kitamura (1969); (5) Robinson 
(1967); (6) Guinan (unpub.); (7) Edalati, Khalesse, & Riazi (1989); (8) This work; (9) Varricatt & Ashok (1999). 



HJD" 




cr\s) Ref. 


HJD" 






Ref. 


HJD" 






Ref 


2410368.9939 


-3114 


610 1 


2436982.9957 


-1608 


240 


1 


2442402.5785 


-1287 


770 


1 


2410562.1139 


-2262 


610 1 


2437378.3104 


-1012 


770 


4 


2442426.4195 


-2477 


770 


1 


2410664.3469 


-2415 


610 1 


2437696.3624 


-2027 


770 


4 


2442426.4225 


-2218 


770 


1 


2411425.4369 


-1648 


610 1 


2437746.3434 


-2065 


770 


4 


2442467.3005 


-3593 


770 


1 


2411993.3909 


-3116 


610 1 


2438089.4114 


-896 


770 


4 


2442785.3675 


-3311 


770 


1 


2412527.3029 


-1451 


770 1 


2438105.3104 


-1258 


610 


4 


2442802.4325 


-1076 


770 


1 


2413242.9557 


-642 


610 1 


2438114.3994 


-1132 


770 


4 


2442802.4345 


-903 


770 


1 


2414333.4539 


-1166 


610 1 


2438384.7384 


-2443 


610 


4 


2442820.5941 


-2240 


240 


1 


2414447.0559 


-491 


400 1 


2438399.5192 


-1272 


520 


4 


2442826.2775 


-1921 


770 


1 


2414878.7180 


-140 


610 1 


2438400.6454 


-2114 


610 


4 


2442826.2865 


-1143 


770 


1 


2415810.2070 


1296 


450 1 


2438406.3387 


-940 


520 


4 


2442835.3485 


-3350 


770 


1 


2416718.9560 


904 


610 2 


2438440.4174 


-902 


770 


4 


2442835.3685 


-1622 


770 


1 


2418309.2921 


2400 


610 1 


2438817.5334 


-2347 


770 


4 


2443161.3775 


-2169 


770 


1 


2419615.6312 


2901 


770 1 


2438818.6687 


-2402 


290 


4 


2443162.5145 


-2077 


770 


1 


2419849.6342 


2811 


450 1 


2438832.3054 


-1936 


770 


4 


2443186.3616 


-2750 


770 


1 


2420138.1572 


2271 


770 1 


2439140.1446 


-2028 


240 


1 


2443202.2725 


-2075 


770 


1 


2420513.0292 


3236 


770 1 


2439164.0002 


-1957 


240 


1 


2443203.3966 


-3106 


770 


1 


2421278.6462 


2557 


450 1 


2439169.6784 


-2087 


770 


1 


2443219.3125 


-1999 


770 


1 


2421648.9832 


4279 


770 1 


2439492.2922 


-1544 


610 


5 


2443430.5776 


-3744 


770 


1 


2422029.5022 


2417 


450 1 


2439518.4104 


-2275 


520 


4 


2443512.3796 


-2519 


770 


1 


2422030.6382 


2422 


400 1 


2439528.6364 


-2057 


770 


4 


2443513.5136 


-2687 


770 


1 


2422558.8492 


2249 


610 1 


2439533.1794 


-2124 


240 


1 


2443587.3586 


-1929 


770 


1 


2422765.5903 


2211 


450 1 


2439802.4034 


-1626 


240 


1 


2443595.2966 


-3103 


770 


1 


2423098.4213 


2212 


610 1 


2439822.8464 


-1967 


770 


1 


2443612.3376 


-2942 


770 


1 


2423406.2533 


1498 


520 1 


2439863.7384 


-2132 


770 


1 


2443880.4246 


-2535 


770 


1 


2423442.6093 


2004 


400 1 


2439870.5314 


-4090 


770 


1 


2443888.3706 


-3018 


770 


1 


2423866.3213 


2492 


770 1 


2439872.8174 


-2870 


770 


1 


2443905.4166 


-2425 


770 


1 


2424667.1393 


671 


610 1 


2439875.0984 


-2082 


80 


1 


2443946.3066 


-2763 


770 


1 


2425052.2353 


1679 


610 1 


2439896.6704 


-3024 


770 


1 


2443971.2946 


-2998 


770 


1 


2425320.3193 


1825 


450 1 


2439904.6324 


-2125 


770 


1 


2444255.2839 


-2668 


240 


1 


2425650.8783 


1816 


770 1 


2439905.7724 


-1774 


770 


1 


2444281.4036 


-3270 


770 


1 


2425990.5203 


1415 


770 1 


2439912.5774 


-2694 


770 


1 


2444606.2986 


-1921 


240 


1 


2426014.3803 


1866 


610 1 


2439912.5847 


-2064 


240 


1 


2444607.4327 


-2080 


240 


1 


2426027.9993 


803 


770 1 


2439912.5924 


-1398 


770 


1 


2444647.1938 


-1810 


240 


1 


2426753.8563 


-53 


770 1 


2439929.6374 


—891 


770 




2444648.3289 


-1883 


240 




2426994.6883 


1009 


610 1 


2439935.3075 


-1730 


240 




2444649.4616 


-2163 


770 




2428596.3576 


242 


770 1 


2439954.5965 


-3632 


770 




2444664.2304 


-2028 


240 




2428922.3748 


403 


770 1 


2440288.5785 


-2330 


770 




2444672.1848 


-1786 


240 




2429301.7763 


133 


240 1 


2440313.5715 


-2133 


770 




2444998.1898 


-2679 


240 




2429308.5903 


-8 


240 1 


2440582.7835 


-2672 


770 




2444999.3295 


-2354 


240 




2429309.7273 


82 


240 1 


2440591.8781 


-2061 


240 




2445015.2389 


-1817 


240 




2429660.7283 


-355 


450 1 


2440964.4665 


-2109 


240 




2445391.2370 


-1704 


240 


7 


2430035.5853 


-687 


520 1 


2440971.2825 


-2079 


240 




2448137.9592"^ 


-449 


80 


8 


2432999.2353 


-2634 


770 1 


2440979.2345 


-2044 


240 




2448608.2433"^ 


-93 


80 


8 


2433367.3203 


804 


770 1 


2440995.1395 


-1887 


240 




2450088.3866 


850 


240 


9 


2434453.2714 


-16 


770 1 


2440996.2715 


-2228 


240 




2450096.3415 


1136 


240 


9 


2434454.4043 


-270 


770 1 


2441725.5335 


-3327 


770 




2450107.6995 


1013 


240 


8 


2434481.6620 


-694 


240 1 


2441765.3075 


-1941 


770 




2450145.1826 


756 


240 


9 


2435515.3604 


-1461 


240 1 


2442059.5105 


-2456 


770 




2450154.2670 


485 


240 


9 


2435534.6759 


-1074 


240 1 


2442092.4525 


-2484 


770 




2450439.3955 


1096 


240 


9 


2436958.0042 


-1675 


240 1 


2442099.2715 


-2194 


770 




2451896.8199 


2036 


150 


8 


2436959.1430 


-1428 


240 1 


2442100.4005 


-2794 


770 




2451945.6648 


1977 


80 


8 


2436977.3169 


-1530 


240 3 


2442116.3045 


-2724 


770 













''Not in the UTC but in the TT scale (see text). 

''The uncertainties have been computed using the procedure outlined in §3.2 
'^From Hipparcos photometry. 
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Fig. 1 . — Fits to the astrometric data (positions and proper motions in the top 4 panels) and LTT curve (bottom panel) for R CMa. Note that the Hipparcos data 
are one-dimensional and thus carmot be represented directly. Instead, we show the Hipparcos position for 10 normal points (epoch groups). Although ground-based 
positions spanning over one centmy have been used to constrain the least-square solution, these are not represented for clarity. For the proper motions, the Guinan 
& laima (1983) estimation and 3 Hipparcos normal points are represented. 
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Fig. 2. — Residuals of the fit to the LTT orbit of the third body. The remaining oscillations have tentatively been modeled with an LTT perturbation caused by a 
fourth body in a circular orbit. See text for fit parameters. 
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Table 3 

ASTROMETRIC AND LIGHT-TRAVEL TIME (LTT) SOLUTIONS FOR THE TRIPLE SYSTEM R CMA. 



Parameter 


Value & Standard error 


TT (mas) 


22.70±0.89 


fia cos 5 (mas yr"' ) 


I68.I±0.7 


fis (masyr"') 


-137.1±1.2 


ai2 (mas) 


1I7.2±5.3 


^12 (°) 


262.9±20.7 


Pl2 (yr) 


92.8±L3 




0.49±0.05 


in n 


91.7±4.7 


0^12 (°) 


10.5±4.3 


r,P™ (HJED) 


2449343±258 


LTT semiamp. (s) 


2574±57 


Mi2 (Mq) 


1.24±0.05'' 


Mi (Mo) 


0.34±0.02 


a, (AU) 


18.7±1.7 




l.I3594197±0.00000010 


T-oEB (HJED) 


2430436.5807±0.0006 



"Adopted from Tomkin 1985. 




-600 -500 -400 -300 -200 -100 
Aa cos8 (mas) 



100 200 



Fig. 3. — Scale projection on the plane of the sky of the orbits of the eclipsing pair of R CMa (small ellipse) and the tertiary component (large ellipse). The 
barycenter of the triple system is marked with a plus sign and the positions of the stars in 2002 are represented as small hollow circles. The inset shows a blowup 
of the region surrounding thfL barycenter. The orbital properties and the sky projection were derived from the simultaneous analysis of eclipse timing residuals and 
Hipparcos astrometry (see §W). 



